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We use ab initio density functional calculations to study the formation and structural as well as 
thermal stability of cellular foam-like carbon nanostructures. These systems with a mixed sp 2 /sp 3 
bonding character may be viewed as bundles of carbon nanotubes fused to a rigid contiguous 3D 
honeycomb structure that can be compressed more easily by reducing the symmetry of the honey- 
combs. The foam may accommodate the same type of defects as graphene, and its surface may be 
be stabilized by terminating caps. We postulate that the foam may form under non-equilibrium 
conditions near grain boundaries of a carbon-saturated metal surface. 
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The last few decades have witnessed an unprece- 
dented interest in carbon nanostructures, the most 
prominent of them being fullerenes[l J, nanotubes [2] 
and graphene [3 . Previously postulated hybrid carbon 
nanofoam structures [4]-[7] with a mixed sp 2 /sp 3 bonding 
character have received much less attention for lack of di- 
rect experimental observation. The growing body of in- 
formation about the formation of carbon nanostructures 
including graphene [8 , nanotubes [9l fTU] and fibers [TT] on 
transition metal surfaces with a particular morphology 
suggests ways that should favor the formation of par- 
ticular carbon allotropes. We propose that previously 
unseen nanostructures including carbon foam may form 
under specific conditions on a metal substrate. 

Inspired by previously postulated carbon foams [4H7] , 
we explore ways to grow such structures on a carbon sat- 
urated metal substrate. We use ab initio density func- 
tional calculations to investigate the equilibrium struc- 
ture, structural and thermal stability and elastic prop- 
erties of the growing system. The foam structures we 
study, which have a mixed sp 2 /sp 3 bonding character 
and resemble a bundle of carbon nanotubes fused to a 
contiguous 3D honeycomb structure, are rather stable 
even as slabs of finite thickness. The foam structure may 
be compressed more easily by reducing the symmetry of 
the honeycombs. It may accommodate the same type 
of defects as graphene at little energy cost, and its sur- 
face may be stabilized by terminating caps. We postulate 
that the foam could form under non-equilibrium condi- 
tions near grain boundaries of a carbon-saturated metal 
surface and should remain stable until T>3, 500 K. 

Our calculations of the equilibrium structure, sta- 
bility, elastic properties and the formation mechanism 
of the carbon foam were performed using ab initio 
density functional theory (DFT) as implemented in 
the SIESTA code[l2]. We used the Ceperley- Alder 
[T3] exchange-correlation functional as parameterized 
by Perdew and Zunger[T4], norm-conserving Troullier- 
Martins pseudopotentials[15 , and a double-^ basis in- 
cluding polarization orbitals. We used periodic boundary 
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FIG. 1. (Color online) Structure and electronic properties of 
cellular carbon nanofoam. (a) Left panels depict individual 
cells of the foam in top and side view. Right panel shows 
the contiguous foam in top view, with individual cells termi- 
nated by different types of caps, (b) Electron density differ- 
ence An(r) in a plane normal to the surface, indicated by 
the dotted line in (a), (c) Side view of the structure of a 
stable minimum-thickness foam slab, (d) Spin density dis- 
tribution p^ — pi in the structure shown in (c), represented 
in the same plane as in (b). The isosurfaces are shown for 
p t - Pi = ±0.05 el./ A 3 . 



conditions for the 3D infinite foam structure and 2D slabs 
of finite thickness. The 3D foam with 10 atoms per unit 
cell was sampled by a fine grid [16] of at least 16x16x16 Ap- 
points in the Brillouin zone. We used a mesh cutoff en- 
ergy of 180 Ry to determine the self-consistent charge 
density, which provided us with a precision in total en- 
ergy of <2 meV/atom. 

The structure of the carbon foam is depicted in 
Fig. [TJa). In top view, it closely resembles the graphene 
honeycomb lattice with two important distinctions. We 
find the optimum lattice constant in the honeycomb 
plane of the foam to be a = 4.81 A, which is about 
twice the graphene value a = 2.46 A. More important, ID 
carbon-carbon bonds in the 2D graphene structure corre- 
spond to 2D walls in the infinite 3D foam structure. The 
foam cells, shown in the left panels of Fig.[]Ja), closely re- 
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semble (6,0) carbon nanotubes. The foam contains 60% 
3- fold coordinated C atoms, labeled sp 2 , and 40% 4- fold 
coordinated C atoms, labeled sp 3 . The gravimetric den- 
sity of the optimized foam structure, p — 2.4 g/cm 3 , 
lies in-between the experimental values [17 for graphite, 
p — 2.27 g/cm 3 , and diamond, p = 3.54 g/cm 3 . We 
find the 3D carbon foam structure to be less stable than 
graphene by AE co h~0A2 eV/atom, which is comparable 
to the Ceo fullerene. 

In Fig. [ljb) we display the electron density difference, 
defined by An(r) = n tot (r) — ^n atorn (r) as the difference 
between the total electron density n to t(?) and the super- 
position of atomic charge densities rw om (r). Charge ac- 
cumulation in the bond region indicates strong covalent 
bonding especially between neighboring sp 2 atoms. Our 
DFT results for the electronic structure indicate that the 
bottom of the conduction band lies below the top of the 
valence band, suggesting that the infinite foam should be 
metallic. In reality, this finding is a well-known artifact 
of DFT that we correct using the LDA+U method, which 
indicates semiconducting behavior of the bulk structure. 

Besides the bulk structure, we also considered and op- 
timized foam slabs of different thickness. We must take 
into account the fact that the surface terminated with 
sp 2 -type atoms, which are shared by two honeycombs, is 
inequivalent to a surface with sp 3 -type atoms, which are 
shared by three honeycombs. The thinnest stable free- 
standing slab, dubbed the 'single-decker' structure and 
shown in Fig. [ljc), has both surfaces of the sp 3 -type. 
It has some commonalities with graphitic nanostructures 
that show magnetic ordering at zigzag edges [T5H22] . Sim- 
ilar to the narrowest zigzag graphene nanoribbon, our 
system displays a flat band near Ep that gives rise to 
spin polarization with antiferromagnetic coupling across 
the slab, as seen in Fig. [TJd). The dominating role of the 
surface reduces the stability of the 'single decker struc- 
ture' by AE co h = 0.95 eV/atom with respect to the bulk 
carbon foam. We note an even-odd alternation in the en- 
ergy as a function of slab thickness in terms of the num- 
ber of hexagon rows, since the slab surfaces may be either 
identical or different. In any case, the role of the surface 
decreases with increasing slab thickness, and reaches a 
much smaller value AE co h~0A6 eV/atom in the 'triple- 
decker', shown in Fig. [ljb) , than in the 'single-decker' 
structure. 

The energy penalty due to unsaturated surfaces may 
be significantly reduced if the slab is attached to a sub- 
strate, or if the cells are covalent ly terminated by caps, 
similar to the dome termination of carbon nanotubes. 
We considered either a hexagon or two adjacent pen- 
tagons as candidate caps to terminate the honeycombs, 
as seen in the right panel of Fig. [TJa). Both caps have 
6 twofold coordinated C atoms at the edge that may 
form covalent bonds with the surface atoms. Assuming 
that all honeycombs on one side are capped and using 
A = 20.04 A 2 for the area of each honeycomb, we esti- 




FIG. 2. (Color online) Defects in the foam, (a) Folding of the 
perfect foam, induced by applying hydrostatic pressure or by 
electron doping. Foam structures containing (b) 5775 and (c) 
558 defects, familiar from defective graphene. 



mated the surface energy reduction associated with cap 
termination to be AE S = —1.03 eV/A 2 for hexagonal 
caps and AE S = —0.25 eV/A 2 for the less-stable two- 
pentagon caps. We need to note that this stabilization 
energy contains the termination energy of both the sur- 
face and the individual unsaturated caps, and that these 
energy terms can not easily be separated. 

Since epitaxy is an issue when considering the pos- 
sibility of foam growth on a metal substrate, we in- 
vestigated the lateral compressibility of the foam struc- 
ture. Our definition is analogous to the elastic response 
of a uniform isotropic 3D structure with volume V to 
hydrostatic pressure P = F/A, given by the force F 
per area A, which is represented by the bulk modulus 
B = — V{dP/ 8V)t- The elastic deformation of the area 
A within a 2D slab structure subject to in-plane hydro- 
static pressure P^d = F/l, given by the force per length 
Z, can be represented by an analogous 2D bulk modulus, 
defined by E>2D = —A(dP2D/dA)T- Of course, we expect 
F>2D to be nearly proportional to the slab thickness. We 
find this value to be quite useful, since it allows to deter- 
mine the critical slab thickness for epitaxial growth on a 
particular incommensurate substrate. 

Applying hydrostatic pressure in the plane of the layer, 
we find that the honeycomb structure may be compressed 
more easily by breaking the honeycomb symmetry than 
by uniformly compressing the honeycombs. The struc- 
ture of the deformed foam, depicted in Fig. |2|a), indi- 
cates the preferential way the foam may fold. For this 
elastic response, we find B^d — 112.9 N/m in the 'single- 
decker' and B2D — 163.9 N/m in the 'triple-decker' struc- 
ture. For the sake of comparison, when considering a very 
thick slab of thickness h, we used the bulk calculation to 
obtain Bsd^B2D /h = 178 GPa. We find this value to 
be much smaller than that of the ideal structure with 
suppressed folding, which had been studied previously [5 
with results similar to our value £?2d/^ = 299.4 GPa. 
Even though the possibility of folding reduces the bulk 
modulus, finite compressibility should still play a signif- 
icant role during foam growth on lattice-mismatched or 
defective substrates. 

Interestingly, we find that foam folding occurs spon- 
taneously when the system is doped by electrons. The 
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structure presented in Fig.[2|a) can be obtained by either 
applying isotropic pressure in 2D or, at zero pressure, by 
doping with 0.2 electrons per C atom. In the latter case, 
we find that folding induced by doping reduced the foam 
energy by 0.19 eV/atom for the bulk structure. 

We also find that the proposed foam structure may 
accommodate a similar type of defects as graphene with 
the main difference that bond rotations in graphene cor- 
respond to wall rotations in the foam. In graphene mono- 
layers, lines of 5775 or Stone-Thrower- Wales [23 , 24j and 
of 558 defects have been observed to accumulate near 
grain boundaries[25-27 and step edges [28 . Their pres- 
ence reduced stress in strained free-standing layers and 
the lattice mismatch energy in adsorbed layers, which in 
this way maintained their epitaxy over large areas. The 
analogous 5775 or 558 defect structures in the foam are 
depicted in Fig. |2|b) and [2|c). Since the foam struc- 
ture is rather flexible, the energy penalty associated with 
these types of defects is relatively small, amounting to 
AE = 0.19 eV/atom for the 5775 structure of Fig. |||b) 
and AE = 0.20 eV/atom for the 558 structure of Fig.|2^c) 
with respect to the perfect infinite honeycomb lattice. 
With a bulk modulus 5^250 GPa, the defective 5775 
and 558 foam structures are slightly more compressible 
than the perfect foam with suppressed cell folding. Sim- 
ilar to supported graphene, these types of defects should 
reduce the lattice mismatch energy on a particular sub- 
strate caused by different lattice constants or, on a poly- 
crystal, across grains with different orientation. 

To find out whether the carbon foam may or may not 
decompose to a more stable allotrope under growth con- 
ditions, we studied its thermal stability by performing 
molecular dynamics (MD) simulations in the tempera- 
ture range 500 K<T<5, 000 K. To avoid artifacts caused 
by small unit cells, we used supercells containing 160 
carbon atoms. For these large unit cells, we used the 
Tersoff bond-order potential [29 in molecular dynamics 
simulations covering time periods of 10 ps using 0.5 fs 
time steps. Our results, presented in the EPAPS on-line 
material [30 , indicate that the infinite foam should be sta- 
ble up to a high melting temperature near 3, 700 K. Even 
though free-standing slabs with finite thickness may be 
thermally less stable, termination by caps or attachment 
to a substrate should increase their thermal stability. 

Inspired by the observed growth of graphene [8] and 
carbon nanotubes[9] on cobalt saturated with carbon, we 
studied possible growth pathways of the foam on this sub- 
strate. To get insight into the foam-substrate interaction 
including optimum lattice registry, we represented the 
Co (0001) surface by a four-layer slab with the two bot- 
tom layers constrained in the bulk geometry. Besides 
the perfect Co (hep) lattice, we also considered fee layer 
stacking when discussing grain boundaries. We consid- 
ered different foam terminations at the interface in or- 
der to find the optimum interface geometry. We found 
that the sp 2 -type terminated foam surface attaches more 




FIG. 3. (Color online) Surface and bulk diffusion of C atoms 
on a carbon saturated Co (0001) surface. Surface diffusion in 
(a) is compared to bulk diffusion in (b) and diffusion along 
a grain boundary in (c). The top panels represent energy 
changes per atom along the optimum diffusion path, which is 
indicated by the dashed line in top and side views, presented 
in the bottom two panels. 

strongly to Co(0001) than the sp 3 -type terminated sur- 
face. The largest reduction of the foam surface energy 
by AE S = -0.75 eV/A 2 occurs, when surface C atoms 
occupy the hollow sites. We should note that this stabi- 
lization energy reflects the reduction of both the metal 
and the foam surface energy. 

Since a realistic representation of the growth mecha- 
nism by molecular dynamics (MD) simulations is cur- 
rently not possible due to time limitations, we discuss 
in the following likely processes that should contribute 
to foam growth and judge their importance according 
to potential energy surfaces. To favor foam growth, we 
need to find a suitable substrate geometry and identify 
growth conditions that promote the formation of foam 
rather than other competing nanostructures[8 , 9 , 31 . As- 
suming that the feedstock are carbon atoms dissolved in 
the substrate, we consider grain boundaries and steps as 
preferential nucleation sites of the foam. Three compet- 
ing processes contribute to the nucleation and growth of 
carbon nanostructures on the surface: surface diffusion of 
carbon, bulk diffusion of carbon inside individual grains, 
and bulk diffusion along grain boundaries that had not 
been considered previously. 

Our results for these three processes are presented in 
Fig. [3l Since surface diffusion of C atoms, depicted in 
Fig. [3^a), does not require displacement of substrate 
metal atoms, it occurs with a low activation barrier of 
only 0.41 eV and should be the fastest process of all. 
The optimum path involves diffusion between the more 
stable hollow sites, with the hep sites being energetically 
favored by 0.28 eV over the fee sites, across less stable 
bridge sites labeled b. 
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FIG. 4. (Color online) Possible formation mechanism of 
the cellular carbon nanofoam, represented by structural snap 
shots in top and side view. Different grains are distinguished 
by color and shading. Initial formation of a graphene nanorib- 
bon along a step edge in (a) is followed by lateral growth of 
honeycomb cells in (b). 

In bulk cobalt, carbon atoms prefer energetically the 
octahedral interstitial sites over the tetrahedral sites. 
The optimum bulk path, presented in Fig. [3jb), involves 
diffusion normal to the surface between octahedral (o) 
sites across barriers at the triangular face centers (fc) of 
the octahedra. We emphasized one triangular face of an 
octahedron by the white dotted line in the middle panel 
of Fig.^b). In this view, the barrier fc site in the center 
of the triangle separates the favored o sites directly below 
and above. Since the Co atoms are closely packed in the 
hep structure, passing through the center of the triangu- 
lar face requires displacing atoms, which requires a high 
activation energy of 3.19 eV. This value is to be consid- 
ered an upper limit, since presence of defects including 
vacancies should reduce the activation barrier for bulk 
diffusion significantly [32]. 

In contrast to a single crystal, the atomic packing at 
grain boundaries is less compact. Consequently, intersti- 
tial carbon atoms may find an energetically less costly 
diffusion path along the grain boundary than in the per- 
fect lattice. A possible grain boundary structure that 
ends in a step edge is shown in the middle and bot- 
tom panel of Fig. [3jc) . The atomic packing in this grain 
boundary resembles that of a simple cubic lattice, with 
interstitial carbon favoring energetically the body center 
be sites in the cube center. The optimum diffusion path 
requires passing through a square face center fc at the 
interface of neighboring cubes. As seen in the top panel 
of Fig. |3jc), the activation barrier for the diffusion along 
this grain boundary is ~1.3 eV, less than half the single 
crystal value. Considering growth conditions similar to 
those in Ref. diffusion to the surface along this grain 
boundary should be ^4xl0 10 times faster than in the 
perfect crystal at T = 900 K according to Arrhenius law. 

With the information at hand about the diffusion rates 
of the carbon feedstock, we proceed to discuss a possible 
growth scenario. The Co structure in Fig. [4] schemat- 
ically depicts three grains, distinguished by color and 
shading. It is plausible to assume that the terrace height 
at both sides of the grain boundary may not be the 
same, yielding a step structure, which is best visible in 



side view. Under growth conditions [8] near 600° C, the 
fastest rate of carbon diffusion to the surface is along 
the grain boundary towards the step edge, where car- 
bon may aggregate to a narrow graphene nanoribbon. 
Since according to our studies a zigzag edge binds more 
strongly to Co than an armchair edge, we consider a 
zigzag graphene nanoribbon attached to the step edge, 
as seen in Fig. j4^a). To best conform to the substrate, 
the nanoribbon acquires a washboard structure, depicted 
in the top panel in Fig.|4^a). The more reactive nanorib- 
bon atoms, which protrude towards the terrace, are more 
likely to form bonds with carbon atoms diffusing along 
the terrace, thus initiating the formation of foam cells. In 
the meantime, atoms or flakes diffusing along the upper 
terrace become the feedstock for the termination of the 
foam layer by caps, as seen in Fig. [4^b). More detailed 
structure information and an animation of the growth 
process is presented in the EPAPS on-line material[30j. 
We hope that this information may encourage follow-up 
experimental studies aiming at synthesizing the carbon 
foam and related carbon allotropes. 

In conclusion, we studied the formation and struc- 
tural as well as thermal stability of cellular foam-like 
carbon nanostructures by performing ab initio density 
functional calculations. We found that these systems 
with a mixed sp 2 /sp 3 bonding character may be com- 
pressed by reducing the symmetry of the honeycomb 
cells. The foam may accommodate the same type of de- 
fects as graphene, and its surface may be be stabilized by 
terminating caps. We postulate that the foam may form 
under non-equilibrium conditions near grain boundaries 
on a carbon-saturated Co surface and should be ther- 
mally stable up to ~3, 700 K. 
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